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Abstract Standing sausage modes are important in interpreting quasi-periodic 
pulsations in the lightcurves of solar flares. Their periods and damping times 
play an important role in seismologically diagnosing key parameters like the 
magnetic field strength in regions where flare energy is released. Usually such 
applications are based on theoretical results neglecting unresolved fine structures 
in magnetized loops. However, the existence of fine structuring is suggested on 
both theoretical and observational grounds. Adopting the framework of cold 
magnetohydrodynamics (MHD), we model coronal loops as magnetized cylin¬ 
ders with a transverse equilibrium density profile comprising a monolithic part 
and a modulation due to fine structuring in the form of concentric shells. The 
equation governing the transverse velocity perturbation is solved with an initial- 
value-problem approach, and the effects of fine structuring on the periods P 
and damping times r of global, leaky, standing sausage modes are examined. 
A parameter study shows that fine structuring, be it periodically or randomly 
distributed, brings changes of only a few percent to P and r when there are more 
than about ten shells. The monolithic part, its steepness in particular, plays a far 
more important role in determining P and r. We conclude that when measured 
values of P and r of sausage modes are used for seismological purposes, it is 
justihed to use theoretical results where the effects due to fine structuring are 
neglected. 
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1. Introduction 


Low-frequency magnetohydrodynamic (MHD) waves and oscillations have been 
abundantly detected in the solar atmosphere with both spectroscopic and imag¬ 
ing instruments since the late 1990s with the advent of the TRACE, SOHO 
and Hinode satellites. This has led to the rapid development of coronal seis¬ 
mology, enabling the inference of atmospheric parameters diff icult t o measure 
directly (for recent revie ws, see e.g., Nakariakov and Verwichte. 12005 : Stepanov, 
Zaitsev, and Nakariakov, [201^. To name but a few, seismological applications 
can offer such key inform ation a s the magnetic held strength in coronal loop s (e.g. ,| 


Naka riakov and Ofman, 200lh and above streamer stalks (Chen et al., 


20111 ). the magni tude of Held-aligned loop flows (Li, Habbal, and Chen, 


2010 


2013 


Chen et aL,[201J), as well a s how the densit y is di stributed transverse to coronal 
loops (Arregui et al., l2007l : Goossens et al.. l200i^ . 

In the context of coronal seismology, sausage waves were examined primar¬ 
ily in connection to their potential for interpreting second-scale quasi-periodic 
pulsat ions ( QPPs) in flare lightcurves (see the review by Nakariakov and Mel¬ 
nikov, l2009|). Two regimes of sausage waves are known to exist, depending on 


the wavenumber k in the direct ion of magnetic structures hosting them {e.g., 
Roberts, Edwin, and Benz,ll98J). When k exceeds some critical kc, the trapped 
regime results, whereby the wave energy is well conhned to magnetic structures. 
If on the contrary, k < kc, the leaky regime results, and sausage waves experience 
apparent damping by emitting waves into the surrounding fluids. For standing 
sausage modes, it is now well established that their period P increases steadily 
with decreasing k until saturating when k is sufficiently small. Identically inhnite 
in the trapped regime for ideal MHD Huids, the damping time r decreases with 
decreasing k unt il a s aturation value is reached for sufficie ntly s mall k {e.g., 
Kopylova et a/., 2007; Nakariakov, Hornsey, and Melnikov, 20121 : Vasheghani 
Farahani et al.. l2014l f . 

Both the period P and damping time t of sausage modes are important 
from the seismological perspective due to their dependence on the atmospheric 
parameters. Let a and L denote the half-width and length of a coronal structure, 
respectively. In addition, let vm {vAe) be the internal (external) Alfven speed. 
For trapped modes, it is well known that the period P is determined by L/vph 
with the phase speed Uph lying i n the r ange between vm and VAe {e.g., Aschwan- 
den, Nakariakov, and Melnikov, 2004[) . On the other hand, for leaky modes the 
period P was found to depend primarily on the internal Alfven transit time 
a/vAi, while the ratio t/P is primarily determined by the density co ntrast pi/pe 
between the structure and its surroundings {e.g., Kopylova et al., 2007l l. For 
both P and r, the detailed form, the steepness in particular, of the transverse 
densi ty distribution can play a subtle role (Nak ariakov, Hornsey, and Melnikov, 
I2OI2I : Hornsey, Nakariakov, and Fludra, l2014l l. Regarding the applications of 
measurements of leaky sausage modes, inferring a/vAi is important given the 
importance for inferring the coronal magnetic Held strength. The inference of 
Pi/Pe and density proHle steepness, on the other hand, is also important because 
these two parameters are key to determining the efficiency of such coronal heating 
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mechanisms as phase-mixing (Heyvaerts and Priest, 1983li and re sonant absorp¬ 
tion (see the review by Goossens, Erdelyi, and Ruderman, 1201 ll . and references 
therein). 

The current conclusions on the behavior of P and r of sausage modes are 
mainly based on theoretical and modeling efforts where the density structuring 
transverse to oscillating loops is either in a piecewise-c onsta nt (tophat) form 
{e.g., Roberts, Edwin, and Benz, Il984 : Kopylova et ai, 200it) or descri bed by 
some smooth functions (e. 5 ., Nakariakov, Hornsey, and Melnikov. l2012h . How¬ 
ever, both theoretical and observational studies indicate that coronal loops may 
not be regularly structured but contain unres olved, fine structuring. For instance, 
densi ty measurements (Tripathi et al, 20091 : Brooks, Warre n, and Ugarte -Urra , 
I 2 OI 2 I I. and temperature diagnostics {e.g., Warren et ai, 20081 : Brooks et aZ.. 12013 1 
suggest that loops with apparent widths > 1000 km are likely to comprise a 
multit ude of fine structures, the scales of which may be down to < 15 km (Peter 
et ai, 201^. On the theoretical side, numerical s imula tions on the temporal 
evolution of coron al loops {e.g., Reale and P eres, ^OOOl : Warren, Winebarger, 


and Mariska. l2003l : Guarrasi, Reale, and Peres. 120101) indicate the need to invoke 
threads with widths ^ 100 km to better reproduce direct observables. Interest¬ 
ingly, further evidence showing the existence of fine struc turing is found from a 
seismological standpoint by Van Doorsselaere et al. (120081 1 when interpreting the 
GoMP measurements of propagating transverse waves. The increasing consensus 
on the existence of fine structuring in coronal structures has also stimulated 
the interest in examining the influence of fine s tructu ring on coronal seismol¬ 
ogy (e. q., M urawski. Nak ariako v, and Pelinovskv. I 2 OOI : Pascoe, Nakariakov, and 
Arber. 12007 : Yuan et al 2015ll. Regarding sausage modes, the study bv Pascoe. 
Nakariakov, and Arber ( 200?! . hereafter PNA071 is of particular relevance since 
it focuses on how fine structuring influences the period P of trapped modes. 
Modeling coronal loops as a magnetic slab with either periodical or random fine 
structuring, PNA07 find that P of fundamental trapped modes is insensitive to 
fine structuring. 

While the study by PNA07 is reassuring for seismological applications based 
on measured values of sausage periods, it is natural to ask: Does this conclusion 
hold also for magnetic cylinders? In addition, what will be the influence of fine 
structuring on the leaky modes, the damping time r in particular? These ques¬ 
tions are addressed in the present study. This manuscript is organized as follows. 
In Section [5] we present the necessary equations and our method of solution. A 
parameter study is then presented in Section [3] to examine in detail the changes 
to P and r due to fine structuring relative to the case where fine structuring is 
absent. Section |4] closes this manuscript with our summary and some concluding 
remarks. 


2. Model Description and Method of Solntion 

We work in the framework of cold (zero-beta) MHD, appropriate for the solar 
corona, and model coronal loops as straight cylinders whose axes coincide with 
the z-axis in a cylindrical coordinate system (r, 9, z). The equilibrium magnetic 
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field B is uniform and also in the z-direction {B = Bz). The equilibrium density 
p is assumed to be a function of r only. In this regard, the modeled fine structures 
are concentric shells sharing the same axis as the cylinder. Considering only fun¬ 
damental standing sausage modes supported by loops of length L, we can express 
the radial component of the perturbed velocity as 6vr{r,z,t) = v(r,t) s\n{kz) 
with the axial wavenumber k being tt/T. It then follows from linea rized , ideal, 
cold MHD equations that (see Nakariakov, Hornsey, and Melnikov, l2012[ l 


dMr,t) 2 , ^ ,15 

dt2 



( 1 ) 


where WA('f) = B/y/4z:p{r) is the Alfven speed. The boundary conditions ap¬ 
propriate for sausage modes are 




Furthermore, the following initial conditions (ICs) 


( 2 ) 


'f'(cOli=o = 


r 

1 -I- ^2 ’ 


dv{r, t) 


dt 


= 0 

t=o 


(3) 


are adopted. While in principle the ICs can be arbitrary, the form given in 
Equation m ensures that only waves with the simplest radial structure are 
primarily excited. 

Similar to PNA07, the equilibrium density p{r) consists of two parts, a back¬ 
ground monolithic one (pmono) modulated by a distribution due to fine structur¬ 
ing (pfs)- In other words. 


p(r) 

= Pmono(?') + PFs{r) , 

(4) 

'inono(^ ) 

= Pe + (Pi - Pe) f{r) , 

(5) 

PFS(?') 

= (Pi - Pe)g{r)f{r) . 

( 6 ) 


The function f(r) decreases smoothly from unity at the cylinder axis (r = 0) to 
zero at infinity, ensuring that Pmono(''’) decreases smoothly from pi at r = 0 to 
Pe at infinity. On the other hand, g{r), which describes the fine structures in the 
form of concentric shells, cannot exceed unity to avoid negative values of p. To 
maximize the effects due to fine structuring, we set the maximum of g{r) to be 
unity. 

With the boundary (Equation ([2])) and initial (Equation ([3])) conditions. 
Equation o can be evolved once /(r) and g{r) are specified. In practice, we 
solve Equation (ED with a finite-difference (ED) scheme second-order accurate in 
both time and space. To save computational time, a non-uniform computational 
grid {ri, * = I, 2, •••,/} is adopted with ri = 0 and rj = lOOOa, where a is the 
cylinder radius. The grid spacing {Ar^ = r^+i — ri,i = I, 2, • • •, J — 1} is iden¬ 
tically 0.001a for r < 4a. Then Ar^ increases in the manner Ar^+i = 1.025Ari 
until Ari reaches yjpijp^ x Ari. From there on Ar^ remains uniform again. 
To ensure numerical stability, a uniform timestep At = 0.8Armin/aA,max is set 
according to the Courant condition, where Ar„iin (^^A,max) denotes the minimal 
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(maximal) value of {Ar^} ({t;a,z = 5/Y^47 rp(r^)}). In response to the initial 
condition (Equation ([3])), disturbances are generated and propagate awa y from 
the cylinder. However, as shown in Nakariakov, Hornsey, and Melnikov (|2012h . 
two regimes can be readily distinguished if one follows the temporal evolution of 
the perturbation at, say, r = a. For k = tt/L larger (smaller) than some critical 
value, this v{a, t) evolves into a harmonic (decaying harmonic) form, correspond¬ 
ing to the well-known trapped (leaky) regime. Numerically fitting v{a,t) with a 
sinusoidal (exponentially decaying sinusoidal) function then yields the period P 
{P together with the damping time t) for trapped (leaky) modes. It suffices to 
note here that P and r depend only on the combination [/(r), g(r); L/a, pi/pe] 
when P and r are measured in units of the internal Alfven transit time a/vAi 
with VAi being Bj^yAirpi. 

A number of measures are taken to ensure the accuracy of the numerical 
results. First, a grid convergence test is made for a considerable fraction of 
the computations, whereby the numerical results do not show any appreciable 
change when the grid spacing is halved. Second, we work only with the part 
of the v{a,t) signal when traveling disturbances have not reached r/, meaning 
that the signal is not contaminated from disturbances reflected off this outer 
boundary. Third, we have tested other forms of the initial condition to make 
sure that the derived values of P and r do not depend on this choice. However, 
since we are interested only in the lowest order modes, we choose not to use a too 
localized v{r, t = 0) since this excites higher-order sausage modes as well. Forth, 
for a number of [/(r), g{r)], we also compute P for trapped modes by formulating 
Equation o into an eigen-value problem and then numerically solving it with 
a M ATLAB boundary-value-problem solver BVPSUITE in its eigen-value mode 
(See Kitzhofer. Koch, and Weinmiilleil 20091 for a description of the solver, and 
Li et al. \ \2Q14 for its recent application to sausage modes). The values of P thus 


found agree remarkably well with those found with our initial-value-problem 
approach (see Figure [2]). 


3. Parameter study 


In this section, we perform a parameter study to see how the fine structuring of a 
multi-shelled loop affects the period P and damping time r of the global sausage 
modes. Given that [/(r), g{r)]L/a, pi/pe] constitutes too large a parameter space 
to exhaust, we choose to fix the density contrast p-Jpe to be 50, in accordance 
with the v alue s uggested by the event reported in Nakariakov, Melnikov, and 
Reznikova ( 2003h . In addition, we fix the form of /(r) to be 


f{r) = exp 



f-Vl 


Va/ . 


p> 1. 


(7) 


Obviously, Pmono(r) is increasingly steep when p increases. This choice of /(r) is 
similar to the symmetric Epstein profile adopted in PNA07 whenp « 2, however, 
a variable p also allows us to see the effects of the steepness of the monolithic 
density profile. Unless otherwise specihed, we fix p at 2. 
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Figure 1. Equilibrium density profiles (p(r)) with periodic fine structuring. The red curves in 
the left (right) column are for p(r) in which fine structuring is given by sinusoidal (cosinusoidal) 
forms. The black curves present the monolithic component, to which a modulation due to fine 
structuring is added. Furthermore, two values for the number of fine structures, 1 (the upper 
row) and 20 (lower), are adopted for illustration purposes. 


3.1. Periodic Fine Structuring 


Let us start with examining fine structuring of the form, 


5(d 


2n7rr 

sin- 

a 

2n7Tr 

cos- 

a 


( 8 ) 


with the integer n representing the number of shells in a loop. Figure [T] shows 
the equilibrium density profile p as a function of r, with sinusoidal (cosinusoidal) 
modulations given by the red curves in the left (right) column. Two values of 
n, 1 (the upper row) and 20 (lower), are adopted for illustration purposes. For 
comparison, the p profile without fine structuring is given by the black solid 
lines. 

Figure [2] shows the dependence on the length-to-radius ratio L/a of the period 
P and damping time r. As labeled in Figure [2 }d, the dashed (dotted) curves 
correspond to the choice of g{r) as a sinusoidal (cosinusoidal) function, and the 
red (green) curves represent the cases where n = 1 (n = 100). We note that with 
the asymptotic Alfven speed v^e defined as i?/v^47rpe, the maximal period that 
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Figure 2. (a) Period P and (b) damping time r of sausage modes supported by loops with 
periodic fine structuring. Here P and r are displayed as functions of loop length L. The black 
curves represent the cases where fine structuring is absent. The dashed (dotted) curves are 
for the results where fine structuring is in a sinusoidal (cosinusoidal) form. Two values for the 
number of fine structures, 1 (the red curves) and 100 (green), are examined. In (a), the blue solid 
line separates trapped (to its left) from leaky (right) waves. Furthermore, the asterisks present 
the periods of trapped modes derived by solving the problem from an eigen-value-problem 
perspective (see text for details). 


trapped modes attain is 2L/vAe- This is given in Figment by the blue solid line, 
which separates the trapped (to its left, where r is identically infinite) from leaky 
(to its right) regimes. For comparison, the black curves represent the cases where 
fine structuring is absent (hardly discernible from the green curves though). The 
curves are found from our fitting procedure, whereas in Figure[2ti- the asterisks on 
the left of the blue line are found by formulating Equation © into an eigen-value 
problem and then solving it with BVPSUITE. Obviously, the periods P found 
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Figure 3. Dependence on the number of shells (n) of 5P and (5r, the fractional changes to 
the periods and damping times of sausage modes supported by thin loops with periodic fine 
structuring. The red (green) curves represent the cases where fine structuring is in a sinusoidal 
(cosinusoidal) form. 



with BVPSUITE agree remarkably well with those from the fitting procedure. 
Figure [5] indicates that while differing in details, the overall tendency for P (t) 
to increase (decrease) with L/a is seen for all the equilibrium density profiles 
considered. In particular, regardless of the profiles, at large values of L/a, both 
P and T saturate at some asymptotic values. In addition, while for n = 1 (the 
red curves) the P and r values show some appreciable difference from the case 
without fine structuring (the black curves), for n = 100 (the green curves) the 
values of P and r can hardly be told apart from the black curves. 

From the computations with the two extreme values of n (n = 1 and 100), 
one may reasonably conjecture that the larger n is, the less the effect associated 
with the periodic fine structuring. To see whether this is true, we employ the 
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fact that neither the period P nor damping time t depends on Lja when L/a 
is sufhciently large. Let Pg and Tg denote the asymptotic values that P and r 
attain, respectively. The n-dependence can then be brought out by examining 
how varying n changes Ps and Tg relative to the cases where fine structuring is 
absent (denoted by the subscript mono). These are defined as 


5P = 


Ps-P 


mono.s 


St = 


Ts - T, 


mono.s 


(9) 


Figure [3] presents (a) SP and (b) 6t as a function of n for g{r) in the sinusoidal 
(cosinusoidal) form as given by the red (green) curves. One can see that indeed 
the fractional variations in P and r tend to decrease in magnitude with increasing 
n. However, even in the cases where n = 1, |i5P| is no larger than 3%, and |i5r| 
remains less than 20%, with the most prominent changes found for |(5 t| when a 
sinusoidal form is adopted for g(r) (the red curve in Figure |3 Jd). This happens 
despite that the equilibrium density profiles with fine structuring are remarkably 
different from the one without it (see Figures [T]a, and[TJ;). The reason for the 
tendency for the periodic fine structuring to have less effect with increasing 
number of shells is related to the transverse profile of the corresponding pertur¬ 
bations. While not shown here, these profiles have spatial scales of the order of 
the cylinder radius. With increasing n, the spatial scale (~ a/n) of the transverse 
density profile decreases {e.g., compare Figures[TJi with[T]D). As one expects that 
the effect of fine structuring maximizes when the two scales are comparable, for 
large values of n, the density fine structuring will have too small a transverse 
spatial scale to influence P or r. 


3.2. Random Fine Structuring 

Compared with the periodic form, a step closer to reality will be a random 
distribution of the fine structures (see also PNA07). This leads us to examine 

9 {r) = ^^, A= ^ A„sin , (10) 

Rmax ^ \ a J 

where A„ and (/)„ are two independent arrays of random numbers in the ranges 
[0, 1] and [0, 27r], respectively. The denominator i?„iax represents the maximum 
that R attains, thereby ensuring that g{r) < 1. Moreover, N stands for the num¬ 
ber of harmonics that take part in the fine structuring, and can roughly represent 
the number of randomly distributed concentric shells. For an illustration of the 
equilibrium density profile associated with such a g{r), see the green curve in 
Figure [Sjo where N is set to be 60. 

Obviously, with increasing N, the relative importance of the harmonics with 
lower n decreases. From Figure[3]we know that for single harmonics, the changes 
to both the period P and damping time r decrease with increasing n. One then 
expects that as N increases, the influence of random fine structuring decreases. 
It turns out that this is indeed the case if one examines Figure HI where the 
influence of fine structuring is represented by SP and Sr, the fractional changes 
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1 10 100 

N 


Figure 4. Dependence on N of 5P and 5t, the fractional changes to the periods and damping 
times of sausage modes supported by thin loops with random fine structuring. Here N repre¬ 
sents the number of harmonics that enter into Equation lllOII , and can be seen as the number of 
randomly distributed shells. At every value of N, each dot represents one of the 100 computed 
realizations of the random structuring. 


to the saturation values Ps and Tg relative to the case where fine structuring is 
absent. Each dot for a given N represents a realization of the random function R, 
and in total 100 realizations are computed. FigureH^ indicates that |5P| is con¬ 
sistently smaller than 6%, meaning that the period of the fundamental sausage 
mode is insensitive to fine structuring. Rather, it is primarily determined by the 
monolithic part of the equilibrium density profile. Note that this conclusion was 
established in PNA07 for trapped modes supported by magnetic slabs, and here 
we have shown that it holds for leaky modes in magnetic cylinders as well. On 
the other hand. Figure shows that the effect due to hne structuring may be 
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Figure 5. Equilibrium density profiles (a) without and (b) with random fine structuring. In 
(a), a series of values for p, the steepness parameter, is shown by curves with different colors. 
In (b), the black curve represents the monolithic component of the density profile given by 
the green curve, which includes contributions from 60 harmonics (see Equation poll'l. For this 
monolithic component, a value of p = 2 is adopted. 


important in determining the damping time: for some realizations St can reach 
47%. However, St becomes less than 10% when N > 10. Let us consider what 
this implies for coronal seismology, assuming that only P and r are available 
as observables and as generally accepted, loops consist of a substantial number 
of fine structures. The negative side is that the information on fine structuring 
can not be probed since the changes to P and r due to fine structuring may be 
well below measurement uncertainties. However, the positive side is that, when 
inverting the observables P and r to infer, say, the internal Alfven transit time, 
it suffices to use theoretical results accounting for only the simpler, monolithic 
equilibrium density profile. 

3.3. Effect of the Monolithic Density Profile 

So far we have fixed the steepness index p to be 2 when specifying f(r). One may 
naturally question what would happen if p is varied? To start, let us illustrate 
how p influences the monolithic part of the transverse density profile as given 
in Figure [S^. As expected, the monolithic density profile becomes steeper as p 
increases. 

Figureshows the variation with p of Pg (the red dots) and Tg (green), the 
saturation values of the period and damping time. A random structuring as given 
by Equation m with iV = 60 is adopted. Any dot at each value of p represents 
one of the 100 realizations of the random function R, while the curves correspond 
to the case where random structuring is absent. The relative changes SP and 
St are then presented in Figures [6 )d and [H:, respectively. One can see that the 
effect of fine structuring tends to increase with increasing p. In fact, for p = 1.1, 
the dots in Figure hardly show any appreciable deviation from the curves, a 
result naturally expected since Figures|6)D and[B]: indicate that |5P| and |(5r| are 
no larger than 1%. In contrast, for p = 100, Figures[Bj3 andjBJ: indicate that |(5P| 
and |(5t| may be up to 3.48% and 7.53%. However, these values at a given p are 
marginal to say the most when compared with the changes to Pg and Tg due to 
variations in the steepness index p. For instance, Figure[B^ shows that in the case 
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Figure 6. Dependence on p, the density profile steepness, of the periods Ps and damping 
times Tg of sausage modes supported by thin loops with random fine structuring. In panel 
(a), the absolute values of Pa and Tg are given, while their fractional changes relative to the 
cases where fine structuring is absent are given in panels (b) and (c), respectively. The random 
fine structuring includes 60 harmonics. For every p value, each dot represents one of the 100 
realizations of the random fine structuring. 
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without fine structuring, increases monotonically with 79 , attaining 4.46 ajvp^^ 
(16 ajvAi) when p = 1.1 (p = 100). In relative terms, with p = 100 is 3.58 
times that in the case where p = 1.1. From this comparison we conclude that at 
least for the density profiles we explored, compared with the fine structures, the 
monolithic part plays a far more important role in determining the period and 
damping time. 


4. Conclusions 


On both theoretical and observational grounds, magnetic loops in the solar 
corona are suggested to comprise a multitude of fine structures with transverse 
scales much shorter than loop widths. From the coronal seismology perspec¬ 
tive, while the effects of fine structuring in the form of multilayered slabs on 
the periods of trapped, standing sausag e mod es were shown to be at most 
marginal (Pascoe, Nakariakov, and Arber, 20071) . it remains to be seen whether 
the conclusion also holds for coronal cylinders and for leaky modes. Working in 
the framework of cold magnetohydrodynamics (MHD), we model coronal loops 
as magnetized cylinders with a transverse equilibrium density profile comprising 
a monolithic part and a modulation due to fine structuring. This kind of fine 
structuring can be thought of as constituting concentric shells sharing the same 
axis with the cylinder itself. We focused on fundamental modes with the sim¬ 
plest transverse structure, namely, the lowest order sausage modes without extra 
nodes between the two ends of the cylinder. The equation governing the trans¬ 
verse velocity perturbation was solved with an initial-value-problem approach, 
thereby the period P and damping time t were derived by numerically fitting 
the signals at some given location away from the cylinder axis. Both P and r 
were shown to saturate at some asymptotic values in the slender-cylinder limit. 
This enabled us to adopt 5P and St, the fractional changes to these saturation 
values relative to the cases where fine structuring is absent, as indicators of the 
effects associated with the presence of fine structuring. 

Starting with an examination of periodical fine structures, we showed that as 
the number of shells increases, the effects on the period and damping time of 
leaky sausage modes decrease. This is attributed to the disparity between the 
characteristic transverse scale of the velocity perturbation and the increasingly 
small scale of the fine structures with increasing shell number. When the shell 
number exceeds ^ 2, both SP and St amount to a few percent. Going a step 
closer to reality, we examined whether random fine structuring can have a more 
prominent effect in determining P and r. While SP and St are indeed larger 
for some realizations of the fine structuring than in the cases with periodic 
structuring, when the number of fine structures (roughly represented by the 
number of harmonics that enter into the random structuring) exceeds ^ 10 , 
neither SP nor St is in excess of 10%. We showed by varying the steepness 
of the monolithic part of the density profile that this change may bring forth 
changes to P and r by a factor of several, far more prominent than the effects 
due to fine structuring. 
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Our findings can be considered positive news for seismological applications of 
the period and dampi ng time of leaky sausage mo des us ing theoretically {e.g.^ 
Kopylova et al., I2nn7l: Vasheghani F arahan i et al. l2014h or numerically {e.g., 
Nakariakov, Hornsey, and Melnikov, l2012h derived P and r where fine struc¬ 
turing is neglected. However, formulating fine structures as concentric shells 
should be considered only a preliminary step closer to reality. Fine structures 
in realistic coronal loops may be organized as randomly distributed strands, 
thereby requiring the problem be formulated as a two-di mensi onal one in the 
plane transverse to coronal cylinders {e.g., Luna et al., 20irJl . Furthermore, 


to simplify our treatment, we have neglected the finite plasma beta and the 
longitudinal variation of the equilibrium magnetic field or densit y. The se two 
factor s are unlikely to be important, though (see Pascoe et al., 120091: Ingiis 
et aZ.. l2009h . 
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